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PACS 87. 16. dp - Transport, including channels, pores, and lateral diffusion 
PACS 41 . 20 . -q - Applied classical electromagnetism 
PACS 05 . 70 . Np - Interface and surface thermodynamics 

Abstract. - We study the non-equilibrium dynamics of lipid membranes with proteins that 
actively pump ions across the membrane. We find that the activity leads to a fluctuating force 
distribution due to electrostatic interactions arising from variation in dielectric constant across the 
membrane. By applying a multipole expansion we find effects on both the tension and bending 
rigidity dominated parts of the membranes fluctuation spectrum. We discuss how our model 
compares with previous studies of force-multipole models. 
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Introduction. — Biomembranes are self-assembled 
structures containing lipids and proteins that form selec- 
tive barriers surrounding cells and organelles. They also 
participate actively in many biological processes, includ- 
ing creation and maintenance of ion gradients through the 
dissipation of free energy pQ. The mechanical properties 
of a fluid lipid-protein membrane in thermal equilibrium 
have successfully been described based on the Helfrich 
Hamiltonian which include quantities such as bending 
rigidity and tension, and for asymmetric membranes also 
spontaneous curvature. From this starting point observed 
properties of both membrane shapes as well as their ther- 
mal fluctuations can be derived [3]. The situation, how- 
ever, becomes different when the membrane is driven out 
of equilibrium by some active process. For these active 
membranes a clear general picture has not yet emerged 
that explains their mechanical properties. 

That non-equilibrium activity can significantly alter me- 
chanical properties has been demonstrated in experiments 
involving ion pumps in lipid membranes. These experi- 
ments include mechanical manipulations by micropipettes 
[2US] as well as observations of shape fluctuations by video 
microscopy [B] . The apparent interpretation of the experi- 
ments is that the activity tends to enhance fluctuations of 
the membrane shape. To explain this enhancement it was 
proposed in [7] that an effect of the ion pump activity is 
to induce a force-multipole in the membrane environment. 
By Newton's 3rd law a monopole is excluded, so a shifted 
force-dipole was chosen as the simplest possibility allowed 



by symmetry. This force-dipole model was further studied 
in [MTT]. 

A shortcoming of jTHTT] is that no specific physical ex- 
planation is given for the origin of the force-dipoles. In this 
letter we propose that these forces can arise due to electro- 
static interactions when an ion is transported across the 
low dielectric membrane interior between the surround- 
ing highly polariscablc water. In the following we calcu- 
late electrostatic forces and consequences for fluctuations 
within a simple model of this process and compare with 
previous work on the force-dipole model. We find that our 
model agrees with a force-multipole picture of the activ- 
ity. However, contrary to the assumptions in the previous 
work [7Hll|. we find that our model leads to a symmet- 
ric coupling between pump diffusion and membrane shape 
dynamics. 

Electrostatic model. — Our model for the mem- 
brane is a surface of thickness 2d with a dielectric per- 
mittivity £2 embedded in a fluid of dielectric permittivity 
e, see fig. [U We choose Cartesian coordinates (a:, y, z) with 
a basis of unit vectors x, y, z. Taking the z-axis along the 
normal direction to the membrane we call the fluid region 
below the membrane, z < —d, region 1. The membrane 
region, — d < z < d, is region 2 and the fluid at z > d 
is region 3. In the membrane region 2 there is a charge 
qo situated at (x c , y c , z c ), which is in the process of being 
pumped across the membrane. We are interested in de- 
riving the forces in our system due to the presence of the 
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charge. To calculate these forces we use the method of 
image charges [12] • By this we obtain the potentials 4>i in- 
each region labelled by i: 
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where a = (e 2 — e)/(e2 + e), & = 2e/(e2 + e)> and rj" and 
7*3 are the distances from the n th image charge in regions 
1 and 3 to the point (x, y, z): 



i n) = ^ {x - x c f + (y - y c y + {z - z^Y 
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with i = 1, 3, = (-l) n jz c - 2nd, and z 3 " ; = {-l) n z c + 
2nd. Finally = rf' is the distance to the actual charge 
situated in the membrane. These potentials give rise to 
forces pointing along the normal direction to the mem- 
brane at three locations: two forces distributed on the 
membrane-fluid interface areas due to discontinuity of the 
Maxwell stresses 

h,2{x,y) = T 2 . zz (x,y, -d) - T 1)ZZ (x,y, -d) (5) 
f2,z(x,y) = T 3tZZ (x,y,d) - T 2 , zz (x,y,d) (6) 

where the Maxwell stresses along the normal direction is 
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and £j = e for i = 1,3. The third force arises from the 
image charges interactions with the real ion: 
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Because of continuity of the lateral components of the 
electric field across the membrane interfaces there are no 
lateral components to the forces. The next step in our 
derivation is to approximate the interface forces fi y2 (x,y) 
and f 2,3(2, y) with point forces at the interface. System- 
atically this approximation corresponds to expanding the 
Fourier transforms of these forces to zeroth order in pow- 
ers of d times wavenumber q. With the convention that 
the Fourier transform of a function g{x, y) is denoted with 
a bar 
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dy e- %{q * x+q ^g{x,y) (9) 
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Fig. 1: The lipid bilayer is portrayed as a band of dielectric of 
permittivity £2 embedded inside a dielectric of permittivity e. 
The dielectric media polarisation, due to the presence of the 
positive ions, give rise to a net negative charge at the interfaces. 
These charges produce forces on the membrane interfaces and 
on the ions. 



and q 



Qx + Qy we obtain 



f 1,2(2, y) w 5(x - x c )S(y - y c )fi,2(z c ) (10) 
f 2,3(2, y) ~ S(x - x c )S(y - y c )f 2 ,z(z c ) (11) 



where 

flA Z c) = flMq=0 
f2,z(z c ) = /2,3|g=0 
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Combining the three contributions for all pumps that are 
in the process of transporting an ion across the membrane 
we find the total force arising due to the activity 
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where (x c>n , y c ,n, z c ,n) are the coordinates of the n th charge 
in the process of going through the membrane. We are 
here assuming a sufficiently low density of pumps such 
that we can neglect interactions between the different ions. 
Eq. (fl"!]) is one of the central results of this paper. It shows 
that our model to a good approximation falls within the 
class of forcc-multipolc models for membrane activity. We 
can therefore use the more general results of [TU] to obtain 
the consequences of the activity for the dynamics of the 
membrane shape. 

Net active force on the membrane surface. In 

the following wc will parametrize the membrane shape by 
the height z = h(x, y, t) of its midplane above the a;y-plane 
at time t. The equation of motion for the membrane shape 
h will be the normal component of the membrane force 
balance equation. Assuming that we are at low Reynolds 
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number with most of the dissipation occurring in the sur- 
rounding water we can write the force balance per area 
as 

/„ +T++T- + / act = (15) 

where / rs is the restoring force associated with elastic 
properties of the membrane, T + and T are the stresses 
from the bulk fluid on each side of the membrane, and 
/ ac t is the additional force induced by the activity. We 
will derive equations of motion for the membrane to first 
order in deviations from a uniform flat membrane. Since 
by symmetry the lateral components of the forces must 
vanish at zeroth order, we are allowed to simply identify 
the normal component of the forces with the z components 
in the following (the difference being of higher order in the 
deviations) . 

The derivation of / ac t starting from a force distribution 
of the form given by eq. (fT4"]) was accomplished pertur- 
batively by a moment expansion in |10j . This derivation 
was achieved by studying the lateral stresses and bending 
moments induced in the membrane by the force distribu- 
tion, and then using that these quantities can be related 
to the membrane force balance (T3"HT5] . The end result 
up to the second moment for this calculation for a nearly 
planar membrane is (10) 



/act • z = a dip Ah + -AQ, 



(16) 



where A = d 2 /dx 2 + d 2 /dy 2 is the 2D Laplacian and 
<7dip and Q are the first (dipole) and second (quadrupole) 
moment of the force distribution 
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Apart from the unitless numerical prefactors then one 
could also have written down eq. ([H>| based on eq. (fT4")) us- 
ing symmetry arguments, linearity and dimensional anal- 
ysis. 

Based on the force distribution (jT4")) we will expand the 
moments as 1161 



Q = n A F( 2) + S 



(19) 
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where is the average concentration of pumps per area, 
nA = nA(x,y,t) is the local concentration difference be- 
tween pumps transporting ions in the positive and nega- 
tive z-direction, and 



F (i) = J dz c X(z c ) -<2/i, 2 (z c ) + df 2 ,3(z c ) + z c f ion (z c ) 

(21) 

F(2)= j dz c X(z c ) d 2 f lt2 {z c ) + d 2 f 2 ^(z c ) + z 2 f ion (z c ) 

(22) 



are averaged moments for a single pump with distribution 
X(z c ) for the ions position in the membrane for a positively 
oriented pump (the negatively oriented having a reflected 
distribution. This reflection is the reason behind the oc- 
currence of the difference «a in eq. (f2TH)). We will assume 
that the overall membrane system is symmetric under re- 
flection of the z-direction such that the average value of 
n A vanishes (up-down symmetry). Finally, S is a zero- 
average noise term. We have neglected noise in Odi p , since 
<7dip enters a term in eq. (|16|) which is already first order 
in deviations from uniformity. Inserting eq. (|19p and (1201) 
in eq. (|16[) we arrive at 



/act • z = n%F {1) Ah + -F (2) An A + ^AS. (23) 

Dynamics and fluctuations. — A closed set of equa- 
tions governing the dynamics of h can now be obtained by 
specifying the remaining forces and a dynamic equation 
for riA- The restoring force / rs can be obtained by func- 
tionally differentiating the free energy F of the membrane 
[17] . Supplementing the Helfrich free energy [5] with cou- 
pling to the riA field we can take F to be 



F = J dA (| {2H) 2 + a + An A H 



X 2 
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(24) 



where H w Ah/ 2 is the mean curvature, dA ps dxdy{\ + 
[(dh/dx) 2 + (dh/dy) 2 ]/2} the area measure, k the bend- 
ing rigidity, oo a tension, A a coupling constant, and the 
parameter x ~ ^bT '/n% for sufficiently dilute pump den- 
sity due to the entropy of mixing. Linear terms in n A 
and h are excluded by up-down symmetry. Functionally 
differentiating with respect to h gives the restoring force 

f rs ■ z = = a Ah - nA 2 h - — An a + 0(h 2 ). (25) 
oh 2 

The motion of the bulk water surrounding a lipid vesi- 
cle will be governed by the low Reynolds number Navier- 
Stokcs equation. From a no-slip boundary condition where 
the water and membrane velocity are matched one can ob- 
tain the stress that the water exert on the membrane. In 
Fourier space this leads to the force [10] 



-Arjqh + Tfrydr 



(26) 



where the bar denotes the Fourier transform, the dot de- 
notes a time derivative, r\ is the bulk viscosity and ?7hydro 
is the thermal noise. The strength of this noise can be 
obtained by applying the fluctuation-dissipation theorem. 

The dynamic equation for the pump density difference 
is a continuity equation. We will write it as 



"A = 



II 



-CIV 



A a Ah + XalT-A 



+ V\\,di«>, (27) 



where Vii is the gradient operator restricted to the mem- 
brane (and thus Vii • Vii = A), O is a phenomenological 
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transport coefficient and V\\,diS ls t ne corresponding ther- 
mal noise. The term in the square brackets of eq. (|27|l 
is the functional derivative of F with respect to n&, ex- 
cept that the index a indicates that A and \ can be mod- 
ified by the activity. That the collective diffusion con- 
stant D a - £lXa can change upon pump activation has 
been observed experimentally for the proton pump bacte- 
riorhodopsin [18J . In the following we will discuss how a 
change of A is related to the active force distribution. 

We assume that the membrane is situated in a quadratic 
frame of area L 2 with periodic boundary conditions. This 
leads to a decoupling of the dynamics of the different 
Fourier modes of h and ri\ . Introducing the column vec- 
tor w = (h(q, t), tia(<7, t)) T , we can write the dynamic 
equations for a single mode in matrix form as 

Bw = -A a w + 77 th + r/a (28) 

where B collects the transport coefficients: 

4r)q 



B = 



o i/ (n q 2 ) 

and A a is a matrix corresponding to elastic coefficients: 



(29) 



A a = 



Kq 4 + a cS q 2 -\q 2 (A - F {2) ) 
-\q 2 K Xa 



(30) 



with c ff = oo + "dip- Note that for a lipid membrane 
the tension ao will adjust when ddip changes to satisfy 
constraints on the total area of the membrane [T9"Il2"0"] . 
The contributions of the thermal noise has been collected 
in the zero average noise vector 77th- The correlations of 
this vector with itself need to be 



m*{q,t)vl{q,t')) = 2k B TBL 2 s(t-t') (31) 



for the fluctuation dissipation theorem to be satisfied in 
equilibrium j^T] . Here 77^ is the transpose as well as com- 
plex conjugate of 77th- The term 77 a contains the noise due 
to the fluctuations of the charge distribution inside the 
membrane 



(32) 



where S is the Fourier transform of the active noise S. We 
will assume that each pump operates independently with 
a typical relaxation time t p . This leads to the correlations 



(S(q,t)S*(q,t'))=L 2 n s T a e 



\t~t 



(33) 



where the noise strength T a will be of the order (.F(2)) 2 
and t^ 1 = t^ 1 + D sp q 2 with D sp being the self-diffusion 
constant of a single pump. 

Let us now return to the issue of the influence of the 
pumping activity on A a . In the absence of activity we will 
have A a = A and A a will be symmetric. This is no coinci- 
dence, since A a is then the Hessian matrix of second-order 



partial derivatives of the free energy F with respect to h(q) 
and n&(q). The equality of the two off-diagonal elements 
of this Hessian matrix is thus an example of a Maxwell 
relation. We will argue that A a is also symmetric in the 
presence of activity, which means that A a = A — F^y 
Our argument is that if we have a non-vanishing F% but 
Sit) — 0, then this situation would correspond to having a 
different static distribution of charges in the active state of 
the pump relative to the passive one. But since there is no 
active noise in the system in this situation, then it would 
just correspond to a thermal system with different con- 
straints on the position of the electric charges. Therefore 
the system can still be described by equilibrium statistical 
physics and consequentially also by a free energy depend- 
ing on the constraints, with the corresponding Hessian 
matrix A a being symmetric. Since A a does not depend on 
dynamic noise then this Maxwell relation must hold also 
when the noise S(t) is turned on. This symmctrization 
of the dynamic equations, with activity modified coupling 
constant A a = A — F^) between pump diffusion and shape 
dynamics, distinguishes the present model from previous 
work on force-multipoles [7H9lfTT]. 

If we solve eq. (|2"8"]l for the equal time correlations of the 
shape h(q, t) (using for instance the same methods as in 
[11]) we obtain 



h(q) 
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where 



k = K-(A-F i2) ) 2 /(4 X a) 
is an effective active bending rigidity, 
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and we have introduced the inverse time-scales 
T D l = D a q 2 , 

T' 1 = (Kq 4 + cr cff q 2 )/(4 m ), 
{kq 4 + cr cffg 2 )/(4 W ). 



(34) 
(35) 

(36) 



(37) 
(38) 
(39) 



In the absence of active noise, r a = 0, then the quadrupole 
moment Fr 2 ) enters the fluctuations only through a mod- 
ification of the bending rigidity. This deviates from the 
expressions given previously in [7H91HT] because the part 
of the dynamical equations corresponding to A a was not 
symmetrized in these papers. In the complete absence of 
activity (Fm\ = i*V 2 ) = T a = 0) we recover the equilibrium 
expression for the fluctuations 
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where 



cq kq 4 + a q 2 
k = k- A 2 /(4 X ) 



(40) 
(41) 



p-4 



Ion pump activity 




Fig. 2: The figure on the left is F( 2 ) (continuous line, in units of 10 -28 Nm 2 ) and F/^ (dashed line, in units of 10 -20 Nm) as a 
function of the position of the charge z°- The figure on the right is «//t e a for A = (continuous line) and A = 1.1 10 -28 Jm 
(dashed line) as a function of the position of the charge z^. The parameters used for both figures are iij = 10 16 m -2 , 
go = e = 1.6 KT 19 C, 2d = 5nm, X = Xa = k B T/n° s , k = Wk B T, e = 80e , e 2 = 2e and T a = (F {2) ) 2 . 



is the effective passive bending rigidity. If we take the 
limit of short wavelengths q — > oo in eq. (|34p assuming 
D sp = D a for simplicity we get 



Hq)\ 
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In the opposite limit q — > we find 
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(42) 



(43) 



Thus the long wavelength behavior is completely con- 
trolled by tension also in the active case, while for the short 
wavelength the active noise (of strength T a ) increases fluc- 
tuations besides the influence of the asymmetry in aver- 
age ion distribution in the pump. This asymmetry en- 
ters through the quadrupolc moment i 7 "( 2 ) of the active 
force, which together with A couples fluctuations in pump 
density to membrane shape fluctuations. Note that the 
dipole moment only influences the shape fluctuations 
through a contribution to the tension cr e fj. In between 
the short and long wavelength regimes there is a highly 
complex behavior influenced by for instance the pumping 
relaxation time r a . 

Estimates of the electrostatic effect. — We can- 
not estimate the magnitude of the electrostatic effects 
by assuming a smeared distribution \(z c ) along the full 
width of the membrane. In this case for instance F^ 
will become infinite due to diverging electrostatic attrac- 
tion when the ion is close to the dielectric discontinu- 
ity at the water interfaces. Instead we will estimate 
the effects by assuming that the ion spends most of its 
time in the membrane around a fixed z c = z° such that 
we have \{z c ) = S(z c — In figure [3] we show the 

dipole moment Fm , the quadrupole moment F^) and 
k/k b s = limq^oo (\h(q) | 2 }/(|Mg)| 2 )eq as a function of z°. 



The plots show that values of k/ n B s ~ 2 — 3 are attain- 
able within our model. These values were obtained for a 
similar quantity called effective temperature (which also 
quantifies the increase in short wavelength fluctuations) 
for the pumps Ca 2+ -ATPase [5] (in this reference they 
also measured A = 1.1 10~ 28 Jm) and bacteriorhodopsin 
[4] (here finding A « 0). The effect at long wavelength 
(small q) is not immediately apparent from the values of 
"dip = n s-^(i) ^ — 10 -4 Nm -1 obtained from figure [H 
The values for Odip are large in magnitude compared to 
the tension change observed by a fluctuation analysis in 
[B] upon pump activation. However, as we have suggested 
in a previous paper |20j . we will claim that a mechanism 
comes into play here that moderates the effect of such in- 
trinsically applied tension: a negative contribution to the 
tension will make the membrane tend to expand. But 
due to the large comprcssional modulus for lipid mem- 
branes they can only expand very little before counter- 
acting elastic forces set in. The situation studied here 
is more complicated than the one studied in [50] though, 
since the tendency of a negative tension contribution to 
increase the area stored in long wavelength fluctuations is 
here competing for membrane area with the increase in 
short wavelength fluctuations (quantified by k/n e g in fig- 
ure 0) . We will study such effects of activity on tension in 
a future publication [53]. 

Conclusions. — In this letter we showed that the elec- 
trostatic interactions of an ion situated within a low dielec- 
tric lipid bilayer fits into the general model of membrane 
activity where the active pumps are modeled as fluctuat- 
ing forcc-multipolcs. This provide a possible microscopic 
picture behind this model. Contrary to the assumptions 
in previous work on the force-multipole model, the mi- 
croscopic electrostatic picture presented here prescribes a 
symmetric coupling between the dynamics of the mem- 
brane shape and the density of ion pumps in the mem- 
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branc. The effect of the mean quadrupole moment of the 
active force is thus simplified so that it acts as a cou- 
pling constant between the protein density and the mem- 
brane curvature. This result might not hold for other mi- 
croscopic mechanisms that generate force multipolcs, e.g. 
steric interactions due to conformational changes of the 
pump during the active process. The electrostatic model 
presented here relies on a number of simplifying assump- 
tions. For instance it ignores interaction with free ions 
in the surrounding water, which must be present for the 
ion pumps to operate. However, the charge of the ions 
inside the membrane is already heavily screened by the 
large jump in dielectric constant at the membrane-water 
interface, so we do not expect this interaction to change 
the strengths of the involved forces by much. For sim- 
plicity we have also assumed a symmetric distribution of 
pumps, i.e., that there arc as many pumps pumping ions 
out of the vesicle as there are pumps pumping inward. As 
a consequence we do not take into account the possibil- 
ity of a significant electric potential building up between 
the two sides of the membrane. Such an electric potential 
would lead to additional effects by modifying the tension 
and bending rigidity }25H27] . Furthermore, the mechani- 
cal properties of the interior of the lipid membrane with 
pumps are assumed here as well as in previous work on 
force-multipoles [7HT0] to be similar to that of an incom- 
pressible fluid. And dielectric properties of the pumps are 
modelled as similar to the remaining lipid membrane. It 
would be interesting to test these last assumptions through 
for instance molecular dynamics simulations of lipid mem- 
branes with ion pumps studying electrostatic interactions 
in the system and the membranes mechanical response to 
such forces. 
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